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CLASSIFICATION  OF  CARTOGRAPHIC  FEATURES 
THROUGH  WALSH  TRANSFORMS 


INTRODUCTION 

In  1980,  a  technique  was  developed  by  the  U.S.  Army  Engineer  Topo¬ 
graphic  Laboratories  (ETL)  for  extracting  a  selected  set  of  cartographic 
features  such  as  road  intersections,  straight-line  roads,  and  rectangular  objects 
using  the  Walsh  transform.  The  technique  was  implemented  separately  in  two 
ways.  First,  a  32-  x  32-element,  solid-state  sensor  array  was  used  to  convert 
optical  energy  of  aerial  imagery  into  an  electronic  signal  that  was  processed 
in  a  minicomputer  to  yield  Walsh  transforms.1  Second,  a  prototype  image 
spectrum  analyzer  (PISA)  was  developed  that  uses  a  large  plasma  discharge 
device  (8.5  x  8.5  inches  lighting  area  with  512  electrodes  each  in  both  x 
and  y  directions)  to  generate  two-dimensional  Walsh  function  patterns 
and  produce  512x512  Walsh  transforms  in  1 4  seconds.2 


The  PISA  produced  successful  results  for  a  selected  set  of  targets 
representing  manmade  cartographic  features  as  stated  above.  Since  the  sensor 
array  system  provides  a  variable  threshold  function,  the  results  were  even 
better.  The  cartographic  features  as  described  with  background  scenes  and 
noise  were  extracted  easily  from  several  aerial  imageries.-1  In  the  above  re¬ 
ferenced  reports,  only  the  extraction  of  manmade  cartographic  features  and 
the  extablishment  of  signal  signatures  for  each  feature  are  described. 


In  this  report,  an  algorithm  was  used  to  classify  and  recognize  a  set  of 
selected  cartographic  features  that  were  decomposed  in  the  Walsh  transform 
domain.  The  method  was  implemented  digitally  with  an  experimental  system 
consisting  of  a  solid-state  sensor  array,  a  minicomputer  as  signal  processor, 
and  a  computer-controlled  translational  state  as  the  image  holder.  Recognition 
examples  for  the  set  of  cartographic  features  with  respect  to  the  window  of 
inspection  are  presented.  Finally,  conclusions  are  given. 


*P.r.  Chen  and  W.W.  Seemuller,  Signal  Signatures  of  Topographic  Features  Using  Analog  Technology, 
U.S.  Army  Engineer  Topographic  Laboratories,  Tort  Belvoir,  Virginia,  ETL-0185,  May  1979.  AD- 
A076  110. 

2P.I-\  Chen,  F.W.  Rohde,  and  W.W.  Seemuller,  Prototype  Image  Spectrum  Analyzer  (PISA)  for  Carto¬ 
graphic  Feature  extraction,  U.S.  Army  engineer  Topographic  Laboratories,  l-ort  Belvoir,  Virginia, 
ETL-0204,  October  1979,  AD-A080  729. 


IMPLEMENTATION 


Digital  Minicomputer  as  a  Feature  Detector  •  An  experimental  system 
was  developed  t hat  will  produce  signal  signatures  < Walsh  Transforms)  from 
the  gray  shade  distribution  of  photo  transparencies  and  that  is  based  on  the 
uniqueness  of  the  signal  signatures  to  detect  and  recognize  the  type  of  carto¬ 
graphic  features  in  a  reasonable  short  time.  The  block  diagram  of  the  system 
is  shown  in  figure  1 . 


A  9-  by  inch  aerial  photo  transparency  is  illuminated  by  a  white 
light  source,  and  a  section  of  the  transparency  is  imaged  onto  a  Reticon  32- 
x  32-element  area  sensor  array.  The  transparency  is  held  by  a  computer- 
controlled,  two-dimensional,  translational  stage  so  that  the  image  can  be 
arbitrarily  positioned  against  the  surface  of  the  array  in  both  x  and  y  di¬ 
rections.  The  array  converts  the  optical  energy  of  the  image  into  a  video  signal. 
The  video  signal  is  digitized  and  then  sliced  in  a  Hewlett-Packard,  2108  mini¬ 
computer  to  become  a  two-valued  binary  signal.  Because  the  value  of  the 
threshold  is  adjustable,  a  very  convenient  means  is  provided  for  isolating 
signals  representing  the  selected  feature  image  from  the  unwanted  background 
noise.  The  "on"  value  of  the  video  is  arbitrary  set  to  “100"  and  ‘‘off" 
to  “0"  and  print,  1  by  a  line  printer  as  a  two-dimensional  binary  array  of 
32  x  32  pixels,  which  represent  the  spatial  signal  signatures  of  the  selected 
topographic  features.  The  sliced  signal  is  masked  with  pregenerated,  two- 
dimensional  Walsh  functions  having  only  two  values,  either  +1  or  -1.  to 
produce  the  Walsh  transform  for  that  given  input  image  or  signal.  This  process 
of  obtaining  two-dimensional  discrete  Walsh  transform  can  be  expressed  as 


,  32  32 

A(i.j)  = 

i  v'  y* 

fix.  y)Wal(i,  x)Wal(j.  y) 

1024  “  ~ 

x=i  y=l 

where  fix.  y)  is 

the  image  binary  array: 

Walti.  x)Wal(j.  y)  is  the  two-dimen 

sional  Walsh  function  of  order  i  and  j;  and  x.  v,  irandj  take  values  from 
1  to  3 2. 4 
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II. I  .  Harniuth.  Sctpuncy  Tht'ttry,  Inundation  and  Application ,  Now  York.  Academic  Press.  Inc.. 
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The  Walsh  transform  coefficients,  also  called  decomposed  signal  signa¬ 
tures,  were  divided  by  a  convenient  constant;  for  this  case.  1,024,  which  is 
the  number  of  pixels  in  a  frame  for  normalization.  The  entire  normalized  set 
of  the  Walsh  transform  coefficients  was  printed  by  the  line  printer.5 


The  detecting  scheme  was  implemented  based  on  the  uniqueness  of  the 
Walsh  transform  of  each  feature  under  consideration.  At  first,  a  reference 
signal  signature  was  established  for  each  cartographic  feature  of  the  entire 
selected  set.  The  test  imagery,  after  being  transformed  into  the  Walsh  domain, 
was  then  compared  to  each  reference  signal  signature  sequentially,  and  classi¬ 
fication  was  made.  Because  cartographic  features  may  appear  in  a  variety  of 
locations  in  the  window  of  inspection  (for  this  case  the  active  surface  of  the 
sensor  array)  and  because  Walsh  transforms  are  neither  translationally  nor 
rotationally  invariant,  two  or  more  reference  signal  signatures  are  required  for 
each  class  of  cartographic  feature  to  avoid  misclassification.  With  this  modi¬ 
fication,  four  out  of  seven  cartographic  features  selected  were  recognized 
without  error,  regardless  of  their  location  with  respect  to  the  window.  The 
rest  ot  the  features  were  also  classified  correctly  in  a  majority  of  the  locations. 
However,  misclassification  occurred  when  these  three  features  were  located 
very  close  to  the  comers  of  the  window.  The  How  chart  for  this  detection 
scheme  is  shown  in  figure  2.  Both  the  magnitude  of  a  single  Walsh  transform 
coefficient  (or  sum  of  a  row,  or  a  column  of  transform  coefficients)  and 
the  ratio  of  each  significant  coefficient  to  A(l.  1)  or  A(2,  1)  or  A(3.  1) 
were  used  as  reference  signal  signatures  for  classification.  At  the  end  of  each 
classification,  the  translational  stage  was  automatically  moved  a  predetermined 
number  ol  steps  in  the  x  and  y  direction,  and  the  new  segment  of  the  test 
image  was  projected  onto  the  surface  of  the  array.  The  same  procedures 
described  above  were  repeated. 
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P  I  .  Chen  and  WAV.  Seeinuller.  Signal  Signatures  of  Topographic  h'carurcs  Using  Analog  Technology. 
U.S.  Army  Engineer  Topographic  Laboralories.  I'ort  Belvoir,  Virginia.  I  TL-018S.  May  1979.  AD- 
A076  110. 
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Detection  Results  •  A  selected  set  of  images  containing  manmade 
cartographic  features,  such  as  road  intersections,  straight-line  roads,  and  rec¬ 
tangular  objects  from  aerial  photo  transparencies  with  a  few  mask  patterns, 
were  used  as  input  test  patterns  for  the  system.  The  spatial  signal  signatures, 
the  associated  Walsh  transforms,  and  the  classification  results  are  shown  for 
each  of  the  test  images.  In  figures  3  through  (>.  the  results  are  shown  for 
a  road  intersection,  a  rectangular  object,  a  horizontal-line  road,  and  a  vertical¬ 
line  road,  all  at  different  positions  in  the  window  of  inspection.  They  were 
detected  and  classified  correctly,  regardless  of  their  location  with  respect  to 
the  window.  In  figure  7,  the  detection  results  are  shown  for  a  diagonally 
oriented  road  intersection.  The  correct  classification  was  obtained  for  a  ma¬ 
jority  of  feature  locations  with  respect  to  the  window.  However,  the  feature 
was  misclassified  when  the  center  of  the  diagonally  oriented  road  intersection 
was  positioned  against  four  corners  of  the  window.  The  misclassification 
occurred  because  the  energy  contributed  by  a  branch  of  the  diagonally  oriented 
road  intersection  is  much  stronger  than  that  of  the  other  branch.  In  figures 
8  and  9.  the  results  are  shown  of  the  classification  for  the  diagonally  oriented 
line  road  having  45  and  135  degrees  angles  with  respect  to  the  x-axis  of  the 
window.  They  are  recognized  correctly  in  a  majority  of  locations.  However, 
the  roads  were  not  recognized  when  they  moved  towards  the  four  corners 
of  the  window.  A  very  dark  area  projected  on  the  array  is  shown  in  figure  10. 
I  or  this  case,  no  recognition  was  possible.  A  few  other  segments  of  the  same 
aerial  photo  transparencies  having  relatively  complicated  scenes  were  also 
tried.  Most  of  them  were  recognized  correctly,  although  some  were  misclassi¬ 
fied  or  not  recognized.  The  results  are  shown  in  figures  11  through  13. 


The  overall  results  indicate  that  the  Walsh  transform  processing  tech¬ 
nique  is  quite  successful  to  detect  well-defined,  linear,  manmade  cartographic 
features.  A  rotational  dimension  may  be  incorporated  in  the  future  to  refine 
the  method  for  detecting  the  above  cartographic  features  at  various  angles  with 
respect  to  the  axes  of  the  window. 
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FIGURE  4.  Continued. 
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MAY  15  19  7?  10  :40 

TRANSFORM 

1  2  3  4  5  €  7  6  9  10  11  12  13  14  15  li 

1  25  2-2  0  0  0  0  1-1-1  0  1  0  0  *  0 

2-16-211  0  5  0  6  0  0  1  0  0  2  0  2  C 

3  -6  0-15-1  6-1-5  0  0  0  1  0  2  0-2  0 

4  -1  0  6  1  -11  0  0  0  1  0-1  0-5  00  0 

5  1  010  0  4  1  0  1-1  0-3  1  0  Of  0 

6  50  -4  0  -1  03  -1  00000  000 

7  -4  0-7  0  -1  -1  -2  -1  l-l  3  0  0  00  0 

8  -1  0  1  0-1  0  0  0  1  0-1  0  0  0  0  0 

$  00  -1  00000  -2  1  000000 

10  0  0  1  0-3  0-1  0  0-1-2  0  0  0  0  0 

11  0000  -2  01020000000 

12  0000400000200000 

13  0040000000002000 

14  10  -2  0000000000010 

15  -  1  0  -3  0  0  0  0  0  0  0  0  0  -1  0  1  -1 

16  00  1  00  000  00  000  0  -1  1 


diagonally  orienteo  road  intersection 

l-’KiURl.  7  Continued. 
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fl*Y  15  197?  10-41 
TRANSFORM 

1  2  3  4  5  €  7  8  9  10  11  12  13  14  15  It 

1  25-fc-l  0  0  0  0-1-2  1  0-1  0  0  0  0 

2  7  -23  -2  0  0  3  0  -2  -1  3  0-1  0  2  0  -1 

3  1  -3  10  5  -4  4  3  0  -1  0  -1  -1  -1  1  1  0 

4  1-3  -4  -7  -6  2  0  0  0  0  2  3-1  0  0  0 

5  0  0  -2  0  10  I  -1  0  0  -1  0  -1  3  1  0  0 

€  0  2  0  0  0  -7  -1  0  2  0  0  -2  0  -2  0  -1 

7  2  -4  3  00  -l  8  1-1  1  1  0-1  02  0 

8  3  -4  0  0  0  -  1  0  -4  -2  0  0  0  -1  0  0  -1 

9  0000-10-2011-11-1000 

10  0  2  0  0  0  2  0  0-2-1  0  1  0  1  0  0 

11  0  0-2-1  1-2  2  0  0  0  1  0  1-1  0  0 

12  0  -1  1  3  2  -2  0  -1  0  0  0  -1  2  -  1  -1  0 

13  00  -1  140000  -1  010000 

14  02  -1  00  -3  0000000000 

15  0010  -1  0100000000  -1 

18  0  -1  00000000000001 

DIAGONALLY  ORIENTED  LINE  ROOD.  45  DECREES 

FIGURE  7.  Continued. 
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Results  for  a  Diagonally  Oriented  Line-Road,  135  Degrees. 
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FIGURE  9.  Continued. 
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FIGURE  10.  Spatial  Signal  Signatures.  Walsh  Transforms,  and  Classification 
Results  for  a  Very  Dark  Area  With  Little  Features. 
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FIGURF  11.  Spatial  Signal  Signatures.  Walsh  Transforms,  and  Classification 
Results  for  a  Horizontal-Line  Road  With  Noisy  Background. 
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FIGURE  12.  Spatial  Signal  Signatures.  Walsh  Transforms,  and  Classification 
Results  fora  Road  Intersection  With  Noisy  Background. 
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FIGURE  13.  Spatial  Signal  Signatures,  Walsh  Transforms,  and  Classification 
Results  for  Small  Streets  With  Noisy  Background. 
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CONCLUSIONS 


1.  The  signal  signature  of  the  spectrally  decomposed  cartographic 
features  is  much  simpler  in  distribution  than  the  spatial  signal  signature  of  the 
same  cartographic  feature  for  all  selected  cases. 

2.  In  most  cases,  the  significant  spectral  components  are  distributed 
among  few,  lower  order  Walsh  transform  coefficients.  Further,  each  transform 
pattern  is  unique  in  itself,  and  it  can  be  easily  distinguished  from  the  rest. 

3.  Two  or  more  reference  signal  signatures  were  required  for  each  class 
of  the  selected  set  of  cartographic  features  since  they  may  appear  in  a  variety 
of  locations  with  respect  to  the  window  of  inspection. 

4.  l  our  classes  out  of  seven  of  the  cartographic  features  selected  were 
detected  and  recognized  without  error,  regardless  of  their  locations  with 
respect  to  the  window.  The  rest  of  the  feature  classes  were  also  classified 
correctly  in  a  majority  of  locations.  However,  misclassifications  occurred  when 
these  features  were  positioned  very  close  to  the  corners  of  the  window.  Nearly 
90  percent  recognition  accuracy  was  obtained  for  the  selected  set  of  the 
cartographic  features. 

5.  A  rotational  means  may  be  incorporated  in  the  feature  to  refine  the 
scheme  to  detect  the  cartographic  features  in  a  variety  of  angles  with  respect 
to  the  axes  of  the  window. 
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